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Abstract 

Aims. We investigate the presence of silicon atoms adsorbed on the surface of interstellar polycyclic aromatic hydrocarbons (PAHs) 
to form SiPAH ;r-complexes. 

Methods. We use quantum chemistry calculations to obtain structural, thermodynamic and mid-lR properties of neutral and cationic 
SiPAH ;r-complexes. 

Results. The binding energy was found to be at least 1.5 eV for [SiPAH]^ complexes whereas it is ~0.5 eV for their neutral coun- 
terparts. From the spectral analysis of the calculated IR spectra, we found that the coordination of silicon to PAH+ does not strongly 
affect the intensities of the PAH^ spectra, but systematically introduces blueshifts of the C-C in-plane and the C-H out-of-plane bands. 
Conclusions. The thermodynamic data calculated for [SiPAH] ^ complexes show that these species are stable and can be easily formed 
by radiative association of Si^ and PAH species that are known to be abundant in photodissociation regions. Their mid-lR fingerprints 
show features induced by the coordination of silicon that could account for (i) the blueshifted position of the 6.2 ^m AlB and (ii) the 
presence of satellite bands observed on the blue side of the 6.2 and 11.2 yum AlBs. From such an assignment, we can deduce that 
typically 1% of the cosmic silicon appears to be attached to PAHs. 
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1. Introduction 

1.1. Astrophysical Context 

Silicon is the most depleted abundant element after iron, repre- 
senting a majo r constituent of interstellar "standar d" dust grains 
(IWhittetll2003l) . essentially in the form of sihcates. ljoiiesl OOOOl) 
concluded from the analysis of depletion patterns and dust pro- 
cesses in different regions of the diffuse interstellar medium 
(ISM) that Si is more easily eroded from dust than any of th e 
other elements. This is in line with the work of lTielensl(ll998h . 
who demonstrated that a significant fraction of the elemental sil- 
icon is locked up in a relatively volatile dust component with a 
binding energy roughly equal to 2.0 eV. 

Polycyclic aromatic hydrocarbons (PAHs) that are amongst 
the dust compo nents offer a large surface for atomic and molecu- 
lar adsorptions. iKlotz et al.l (1995) have established a correlation 
between metal-PAH binding energies and the observed depletion 
of these atoms in the ISM. This model predicts that the abun- 
dance of complexed PAH molecules can explain 5-10% of the 
depletion of metallic atoms in the gas-phase of the diffuse inter- 
stellar medium. This work also emphasizes that Fe and Si are the 
most abundant metals bonded on PAHs, if Si is considered as a 
metal although it is not in strict chemical terms. 

Interstellar PAHs are revealed by their emission in the aro- 
matic infrared bands (AIBs) which set aglow the general inter- 
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stellar medium (iLeger & Pugetll 19841: lAUamandoIa et al.|[T985l) . 
The interpretation of the AIBs has motivated a l ot of observa- 
tional, experimental, and theoretical studies (see lTielensll2008l 
for a recent review), but it is not yet entirely resolved, leaving 
the door open to search for new types of PAH-related species 
such as SiPAH 7r-complexes. 



1.2. Silicon-PAH in Physical Chemistry: Experiments and 
Calculations 



ISrinivas et all (Il992h and iJaeger et alJ (l2005h demonstrated 
theoretically and experimentally that the ;r-system com- 
plex [SiCeHg]^ is more stable than all the other isomers. 
Experimental studies on [SiPAH]^ complexes are scarce. In the 
selected ion flow tube (SIFT) experiments caiTied in Bohme's 
group, it was clearly established that the capture of ground state 
atomic silicon ions Si^ (^P) by benzene and naphth alene is an 
efficient process in the gas phase dBohme et alJll99i]) . These au- 
thors also suggested a possible chemical role of the large PAH 
molecular surface in the catalysis of small S i-bearing molecules 
of astrophysical interest dBohme et al.lll989i) . The radiative as- 
sociation of Si^ with a variety of PAHs was studied in the 
low-pressure conditions of a Fourier transform ion cyclotron 
r esonance (FT ICR) mass spectrometer in the group of Dunbar 
(IDunbaret aI.lll994t jPozniac & Dunbaiil l997i) . Charge transfer 
and condensation reactions (insertion of Si between C and H 
with loss of H) were found to compete efficiently with the asso- 
ciation process. The results obtained in these low-pressure con- 
ditions can be a good indication of the processes that are likely 
to occur in the ISM, that is to say that [SiPAH] ^ complexes can 
be formed efficiently and play a role in the chemistry there. 
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Our work aims at investigating the possible relevance of 
[SiPAH]^ complexes in the ISM, emphasizing two aspects: the 
stability of these systems and the influence of the coordination 
of Si on the IR spectra of PAHs. We report extensive calcula- 
tions of structural, thermodynamic, and mid-infrared spectro- 
scopic properties of ground state [SiPAH]^ complexes for four 
compact PAHs of increasing size : naphthalene CioHg, pyrene 
CieHio, coronene C24H12 and ovalene C32H14. Using the meth- 
ods described in section 2, we present the Si-PAH^ binding en- 
ergies and structures in section 3.1, and the IR calculated spectra 
in section 3.2 where systematic effects on band intensities and 
positions are also discussed. In section 4, the astrophysical im- 
plications are discussed and our main results are summarized in 
section 5. 



2. Methods 

Geometry optimizations and frequency calculations were per- 
formed using well-established quantum chemical techniques 
in the framework of D ensity Functional Theory (DFT, 
[Jones & GunnarssonI 11989! ) with the Gaussian03 package 
tFrisch et al. 2003i) . 

The hybrid B3LYP functional was chosen according to its 
wide use in the study of the ground-state properties of PAHs and 
related molecules (see for in stance Langhoff 1996; Bauschlicher 
l2002l:lBauschUcher et all2008l) . The choice of the D95 basis set 
augmented with a complete set of diffuse and polarization func- 
tions, hereafter indicated by D95-H-**, was made after a com- 
parative study with other basis sets and different combinations of 
diffuse and polarization functions. Geometry optimizations were 
systematically followed by the full vibrational analyses to con- 
firm that the geometries are minima on the potential energy sur- 
face. These were also used to determine the discrete IR spectra 
in the harmonic approximation. 

Errors in the calculated values of the binding energies can 
come from uncertainties inherent to DFT methods and the finite 
character of the basis set. To give some insights into the rele- 
vance of the DFT choice and especially of the functional and the 
basis set choices to calculate binding energies, some test calcu- 
lations were performed on the ionization potentials (IPs) of sili- 
con and PAHs (cf Tab. 1). The calculated values of the IPs were 
found to be lower than the experimental ones. For PAHs, the av- 
erage of this discrepancy is about -0.26 eV (from -0.20 to -0.33 
eV) for vertical IPs. This level of calculation provides the right 
relative order compared to experimental data, with IPs system- 
atically underestimated by 4-5%, which is sufficiently accurate 
for the purpose of this work. 

In the study of the vibrational properties of [SiPAH]^ com- 
plexes, scaling factors have to be applied to the calculated har- 
monic frequencies at the DFT level for basis set incompleteness 
and anharmonicity effects. We determined the scaling factors 
for PAH*^ by comparing the calculated spectrum of each PAH° 
to the experimen tal spectrum recorded in co ld matrix of argon 
for naphthalene dHudgins & Sandfordlll998h and i n cryogenic 
matri x of neon for pyrene, coronene and ovalene dJoblin et all 
I1994I) . We subsequently applied a scaling factor derived from 
the average of the scaling factors determined for each PAH to the 
calculated spectra of PAH^ and the related [SiPAH]^ complexes 
in each spectral region of interest. Three scaling factors were 
used, one for the 10-15 fim (1000-700 cm"') range, one for the 
5-10 yum (2000-1000 cm"') range, and one for the 3 fj.m (3300 
cm"') C-H stretch range where modes of larger anharmonicity 
are present. 



Table 1. Computational and experimental ionization potentials 
(IPs) of the four considered PAHs and of Si. Computational val- 
ues are obtained at the B3LYP/D95-I--I-** level of theory. 





IP (eV) 
calc." 


exp.* 


AlF (eV) 
calc. 


exp. 


naphthalene 


7.82/7.91 


8.14 


-0.30/-0.21 


-0.01 


CioHg 










pyrene 


7.10/7.16 


7.43 


-1.02/-0.96 


-0.72 












coronene 


7.02/7.09 


7.29 


-1.10/-1.03 


-0.86 


C24H12 










ovalene 


6.40/6.43 


6.71 


-1.72/-1.69 


-1.44 


C32H14 










Si 


8.12 


8.15 







adiabatic/vertical IF. * Experimental values are extracted from 



the NIST database. ' MP (eV) = //'(PAH) - //'(Si). 



3. Results 

3.1. Thermodynamics and Structures 

We present in this part thermodynamic and structural properties 
of [SiPAH]^ complexes, completed with thermodynamic prop- 
erties of model [SiPAH]" complexes. The calculated relative 
enthalpies at K (Ai/(0 K)) of the lowest energy isomers of 
[SiPAH]° and [SiPAH]+, with respect to the most stable dissoci- 
ation products (Si(^P) + PAH) and (Si(^P) + PAH+) are reported 
in Tab. 2. Binding energies (Eh) correspond to -A//(0 K) values. 
A structure is considered stable when its E^, is positive. In the 
case of the cationic complexes, the relative enthalpies of the re- 
action Si+(2p) -1- PAH ^ [SiPAH]+ at K (AH'{0 K)) are also 
reported. 

All [SiPAH]"^^ complexes were found to have positives Eh, 
providing evidence for their thermodynamic stability. For the 
[SiPAH]^ complexes, the Eb of the most stable structures were 
found to lie between 0.49 and 0.58 eV, leading to an average 
value of ~0.5 eV. They were found to be in a doublet spin-state. 
For the [SiPAH]^ complexes, the Et of the most stable struc- 
tures, found to be in a triplet spin-state, decreases from 2.67 eV 
to 1.69 eV as the size of the PAH increases. Interestingly, Eh de- 
creases when AIP increases. This is consistent with the covalent 
nature of a bond, assuming the same overlap between the two or- 
bitals in interaction. When AIP is equal to zero. Eh is maximal. 
This is the case for [SiCioH8]+ (A/P(exp) = -0.01 eV). As the IP 
of PAH decreases roughly exponentially with the size, approach- 
ing a value of 6.0 eV, the value of Eh is expected to decrease to 
a value of roughly 1.5 eV for larger species. An examination of 
the Mulliken charges on the silicon atom and the PAH molecule 
shows for the larger species that the positive charge is mainly 
located on the PAH. This is consistent with the AIP values. 

In the case of the formation of a [SiPAH]^ molecule from a 
Si^ and a PAH, we determined that that is an exothermic process 
liberating more than 3.0 eV on average. 

All possible ;7r-interaction sites of Si with PAHs were con- 
sidered to find the lowest-energy isomers. Their structures are 
displayed in Fig. 1. We found two low-energy isomers for 
[SiCieHio]^ and three low-energy isomers for [SiC32Hi4]^ with 
very close binding energies: 2.02 and 1.86 eV for the two iso- 
mers of [SiCi6Hio]+ and 1.69, 1.58, and 1.57 eV for the three 
isomers of [SiC32Hi4]^. 

In all the most stable optimized structures Si preferentially 
coordinates on the edge of the PAH, mostly in a two clos- 
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C.[SiCi6Hior(2) D. [SiC24Hi2]^ 




Figure 1. Lowest energy stable structures for the cationic complexes of silicon with naphthalene CioHg (structure A), pyrene CieHio 
(structures B, C), coronene C24H12 (structure D), and ovalene C32H12 (structures E, F, G), optimized at the B3LYP/D95++** level 
of theory. The values expressed in A indicated on each structure correspond to the Si-C bond lengths. 
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Table 2. Electronic states and calculated relative enthalpies at K of the lowest energy isomers of [SiPAH]"^^, with respect to the 
most stable dissociation products (Si(-'P) + PAH) and (Si(-'P) + PAH+) (corresponding to AH(0 K)) and with respect to (Si^(^P) + 
PAH) (corresponding to AH'(Q K)) at the B3LYP/D95++** level of theory. £i(Si-PAH"/+) = -AH{Q K) 





electronic 
state 


AH(0 K) 
(eV) 




electronic 
state 


AH(0 K) 
(eV) 


AH'(0 K) 
(eV) 


[SiCioHs]" 




-0.54 


[SiCioHg]^ 


2A' 


-2.67 


-2.99 


[SiCieHio]" 


3A" 


-0.58 


[SiCi6H,o]+ (1) 


2A 


-2.08 


-3.11 








[SiCi6H,o]+ (2) 


'A 


-1.86 


-2.89 


[SiC24Hi2]° 


3A" 


-0.58 


[SiC24H,2]+ 


2A 


-1.95 


-3.09 


[SiC32Hl4]° 


3A" 


-0.49 


[SiC32Hl4]+ (1) 


2A 


-1.69 


-3.40 








[SiC32H,4]+ (2) 


'A 


-1.58 


-3.29 








[SiC32H,4]+ (3) 


2A 


-1.57 


-3.29 



est neighbour C-atom configuration for neutral species and in 
a three closest neighbour C-atom configuration for cationic 
species (the so-called 772 and 773 edge coordination, respec- 
tively). Interestingly, the 7/2 coordin ation mode was also f ound 
in the case of [FePAH]" complexes (I Simon & JobUnll2007l) and 
[CuPAH]^ complexes (Dunbar 2002), and the 773 coordinatio n 
mode in the case of [SiCeHs]^ complex (Jaeger et al.l l2005l) . 
These coordination modes imply both a curvature of the edge 
of the C-skeleton and ring distortions which had been akeady 
observed in the case of [TM-CfiHsl'^ (TM: transition metals) cal- 
culated structures (Jae ger et al.l2004l) . 

Besides the fact that silicon is mainly present in its ion- 
ized form in the gas-phase of photodissociation regions (PDRs) 
where PAH"^^ are largely detected, the cationic complexed 
species are more interesting to study for their higher stability. 
Therefore only the spectral features of the latter species and 
their astrophysical implications are discussed in the following 
sections. 

3.2. Influence of the Coordination of Silicon on the IR 
Spectra of PAH+ 

As mentioned in Sect. 1, much work has been devoted to the 
identification of the caiTiers of the astronomical AIBs. The main 
bands are located at 3.3, 6.2, 7.7, 8.6, 11.2 and 12.7 fim. The 
3.3 fim band is attributed to the C-H stretch (vc-h), the 6.2 yum 
band to pure C-C stretches (vc-c), the 7.7 fim band to a mixture 
of C-C stretches and in-plane C-H bends, the 8.6 fim band to 
pure in-plane C-H bends (6c-h)- The 1 1.2 and 12.7 fj.m are both 
assigned to out-of-plane (oop) C-H bends (jc-h)- The positions 
of these yc-H reflect the number of adjacent C-H groups on the 
edge of the PAH C-skeleton. A classification was established by 
iHonv et alj dlOOl). The C-H groups carried by aromatic rings 
with no other C-H group are called solo C-H groups and lead to 
IR activity around the 1 1.2 fim AIB. Likewise, C-H groups with 
one and two adjacent C-H groups are respectively called duo and 
trio C-H groups and lead to IR activity in the range of the 12.7 
fj.m AIB. 

The IR calculated spectra of all stable isomers of [SiPAH]^ 
complexes are reported in Fig. 2 in the [3-15] fim range. The 
spectrum of the corresponding bare PAH^ is reported at the bot- 
tom of each figure to make comparisons easier. After calibration 
the calculated discrete spectra are convoluted by a Lorentzian 
line profile with an assumed FWHM (full width at half maxi- 
mum) of 0.2 fim, which accounts reasonably well for the band 
broadening of the astronomical spectra. The h-Wfim/ho-Hiim in- 
tensity ratios which are often used to analyze observed spectra 



are presented in Table 3, along with the relevant shifts occuring 
for the jc-H and vc-c bands. We emphasize that the extracted 
shifts induced by the coordination of silicon to the PAH depend 
only marginally on the spectral calibration (a maximal deviation 
of -0.02 fim without applying scaling factors). 

As can be seen in Fig. 2, the global features of the IR spec- 
tra of [SiPAH]^ complexes remain close to those of bare PAH^. 
This contrasts with the effect of the coordination of Fe to PAH^ 
(Hudgins et al. 2005; S zczepanski et alJ 120061: ISimon & JobUnI 
2007: Simon et al.t2 008) and can be understood considering that 
the charge remains mainly located on the PAH with respect to 
the values of the IPs of Si and PAH. In particular, the 3.3 fim 
band is hardly affected by the coordination of silicon, in either 
position or intensity. 

For the small [SiCioHg]^ complex, the h-miim/ho-iS/jm ra- 
tio decreases substantially with the coordination of silicon (2.92 
to 0.92) whereas for the largest species, the h-Wfiml ho-i5nm ra- 
tios remain of the same order of magnitude (2.86 to 2.57 for 
pyrene, 5.24 to 2.87 for coronene, 7.27 to 5.66 for ovalene; tak- 
ing the average value for several isomers). A simple explanation 
in terms of partial charge transfer seems sufficient to account for 
this trend. Given the fact that PAH^ molecules s how strong fea - 
tures in the [6-10] fim range while PAH" do not (Langhoff 1996), 
the observed decrease of the h-ioiim/ho-is/jm ratio with respect 
to the one of bare PAH^ is due to the coordination of silicon 
which reduces the charge carried by the PAH molecule. When 
the size of the PAH grows, the charge remains mainly located on 
the PAH with respect to the values of their IPs, therefore the co- 
ordination of silicon has a moderated effect on the total IR activ- 
ity in this region, in particular for the largest species. This differs 
with the effect of the coordination of Fe on PAH^ which leads to 
a stron g decre ase of the h-Wfim/ho-iSiim ratio (Simon & Jobli^ 
120071: ISimon et arjr2008i) . 

In contrast, the coordination of silicon leads to significant 
wavelength shifts. Several new modes are activated, mostly in- 
tense in the [6-8] fim range, while the activity in the [8-10] fim 
is reduced. All the coiTesponding vibrations in this region are 
initially in-plane features before complexation. Since the coor- 
dination of silicon on PAH surface breaks the plane of the PAH 
molecule, that leads to a scattering of the molecule vibrations in 
many coupled modes. We point out the appearance of strong IR 
activities in the 6.0 fim range, due to vc-c modes with C atoms 
bounded to Si. Then [SiPAH]^ show some intense vc-c bands 
while introducing an IR activity at shorter wavelengths than the 
band of the bare PAH+ with an apparent blueshift average value 
of -0.17 fim. However, the value of the shift seems to decrease 
when the size of the complexed PAH molecule increases, consis- 
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Figure 2. IR spectra of the most stable structures of PAH^ and their related Si-complexes calculated at the B3LYP/D95-I--I-** level 
of theory. Discrete spectra are convoluted by a Lorentzian profile with a FWHM of 0.2 fim. The derived average scaling factors are 
0.953 for the 3-5 /zm range, 0.976 for the 5-10 yum range, and 0.966 for the 10-15 yum range. 
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Table 3. h-ioiiml Im-xSfim ratios (/?) for the ground states of PAH^ and their complexes with Si and the effects of the complexation 
on this ratio. 





p 
K 


/v^alr7\rl )/ 
i\yrJ\ri ) 


A oop C-ri bend 


OA oop t^-rl Dcnu 
olr7\rl VS i Al 1 


A ip L.-v^ SLreicn 


OA ip siicLcn 

Olrrtrl VS 1 A 1 1 




2.92 




13.28 




6.60 (6.55) 




[SiCioHs]^ 


0.97 


0.33 


13.19 


-0.09 


6.43 


-0.17 








12.02 


new 








2.86 




11.46 




6.47 




[SiCi6H,o]+ (1) 


3.56 


1.24 


11.68 


+0.22 


6.22 


-0.25 








11.15 


-0.31 ; new 






[SiCi6Hio]+ (2) 


1.58 


0.55 


11.85 


+0.37 


6.18 


-0.29 


C24Hl2^ 


5.24 




11.41 




6.38 




[SiC24H,2]^ 


2.87 


0.55 


11.48 


+0.07 


6.21 


-0.17 








11.32 


-0.09 ; new 






C32H14* 


7.27 




11.11 




6.23 




[SiC32Hi4]+ (1) 


5.01 


0.69 


11.03 


-0.08 


6.13 


-0.10 


[SiC32Hl4]+ (2) 


6.63 


0.91 


11.49 


+0.38 


6.11 


-0.12 








10.79 


-0.33 ; new 






[SiC32Hl4]+ (3) 


5.34 


0.73 


11.15 


+0.04 


6.13 


-0.10 








10.95 


-0.16 ; new 







* All values are expressed in iim. 



tently with the reason that induces this shift: the larger the PAH 
is, the less significant the break of its plane is. The calculated 
value of -0.17 /im has therefore to be taken with caution. 

In the [10-15] yum range of jc-h modes, the coordination of 
silicon leads to a splitting of the bands, but the total IR activity is 
preserved. Even though the jc-h band positions depend on the 
number of adjacent hydrogens which varies from four in naph- 
thalene to one in ovalene, complex formation increases the IR 
activity in the shortest wavelength domain. Indeed, the coordina- 
tion of silicon causes a splitting of the main jc-h band of PAH^ 
in two bands: one at lower, and the other one at larger wave- 
length with respect to those of PAH^. For instance, the bands 
at lower and larger wavelengths shift respectively towards 1 1 .43 
and 11. 97 jjm for [SiCifeHio]^ (1), while the corresponding band 
of CieHio^ is located at 11.76 jim. The band at lower wave- 
lengths, located in the solo domain, is due to the jc-h mode 
of the C-H bond with its motion hindered by the presence of 
Si. An interesting feature is that this band is intense and occurs 
for all [SiPAH]^ complexes with Si coordinated above a C atom 
bounded to a H atom (cf Fig. 1, [SiCioH8]+, [SiCi6Hio]+ (1), 
[SiC24Hi2r, [SiC32Hi4]+ (2) & (3)). The average shift of this 
band position with respect to the position of the main jc-h band 
of the corresponding bare PAH^ is -0.22 /im. Contrary to the 
shift observed in the 6.0 //m range, this phenomenon of splitting 
does not depend on the molecular size but only on the geometry: 
the coordination of Si above a C atom bounded to a H systemat- 
ically hinders the C-H vibration. For all these isomers, this band 
is the most intense when its intensit y is normaliz ed to the num- 
ber of involved H atoms (cf Fig. 3). IJoblin et al.l (119941) pointed 
out from experimental spectra on neutral PAHs the four time 
stronger intensity of solo yc-H than other jc-h, which is in line 
with our calculations. The strong intensity and the position of the 
new jc-H induced by the coordination of silicon are consistent 
with a solo character. 



4. Astrophysical Implications 

Our thermodynamic results show that the binding energies of 
[SiPAH]^ complexes are at least 1.5 eV. Since they are three 
times larger than those of [SiPAH]*', [SiPAH]^ complexes are 
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Figure 3. Comparison between the jc-h intensity per H atom 
for PAH^ (rings) and [SiPAH]^ (diamonds for jc-h additional 
band induced by the coordination of Si and crosses for the re- 
maining yc-H, strongest band in case of isomers) considered as 
a function of the system size Np. The square represents the solo 
jc-H intensity per H atom for C32Hi4^. 



more likely to survive in the conditions of the ISM than their 
neutral counterparts. Furthermore the formation of [SiPAH]^ 
complexes by radiative association is an exothermic process of 
3.0 eV without any expected activation barrier (cf Tab. 2, cal- 
culated Ai/'(0 K)). The formation of [SiPAH]+ could therefore 
occur in many regions where Si^ and PAH*^ are abundant, as in 
PDRs (see for instance Kaufman et al. 2006). These [SiPAH]+ 
complexes could then represent a significant fraction of PAH- 
related molecules in such regions. 
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Figure 4. The IR spectrum of the Red Rectangle nebula retrieved from the ISO archive ; (a) the [5.8-6.8] fim range and (b) the 
[10. 8-11. 8] fim range with their underlying continua (dotted line). 



The spectral analysis of the calculated spectra of [SiPAH]^ 
provides insights for astronomical detection of such species. 
Indeed, the coordination of silicon to the R\H^ surface has two 
main effects: 

1. a strong IR activity at the blueshifted side of the bare 
PAH^ vc-c band. The corresponding shift is going from -0.30 
down to -0.10 fim, and seems to be decreasing when the size of 
the PAH increases. 

2. the occurrence of an additional jc-h band for many iso- 
mers, blueshifted with respect to the lowest wavelength jc-h 
band of PAH^. The average shift is -0.2 ftm, independent on the 
size of the complex but characteristic of the coordination geom- 
etry. 

Following th e profile variability of the 6.2 fim AIB from 
source to source. [Peeters et aTl ( |2002|) pointed up two main com- 
ponents at 6.2 and 6.3 fim - classes A and B. This has moti- 
vated further expe rimental and theoretical studies. In particular 
iHudgins et alj (1200 5) computationally established that the class 
B can be explained by the emission of large PAH^, while the as- 
signment of the class A feature is not that clear: the effect of the 
substitution of one or more C atoms by hetero-atoms (N, O, Si) 
and the coordination of metal ion (Fe^, Mg^, Mg^^) to PAH was 
investigated. That study puts forward the role of the insertion of 
N atoms wi thin the C- s keleto n to resolve this question. On the 
other hand, iPino et al] ( l2008h recently showed on the basis of 
IR spectroscopy of soot-like material that the position of the 6.2 
fim AIB can be explained by an evolutionary scenario of car- 
bon dust through its aromatic to aliphatic content. This proposal 
awaits the right analogues in the laboratory of the "astro-PAHs", 
which have to be small enough to respect the stochastic heating 
constraint. 

The vc-c shift induced by the coordination of Si to PAH^ 
that we calculated is at least -0. 1 fim. For larger PAHs, we ex- 
pect a similar shift due to the preference of Si to coordinate on 
the edge of the PAH, whatever its size, and then its abilility to 
break its plane. Therefore this value can be applied to interstel- 
lar PAHs which are l ikely larger than the spe cies we studied 
jLeger & Pugetlll984tlBauschlicher et alj|2008l) . Our study pro- 
vides an alternative scenario to the polycyclic aromatic nitrogen 
heterocycles (PANHs) explanation and the aliphatic to aromatic 
evolution scheme for the shortest wavelength position of 6.2 fim 
AIB. 



Considering that the 11.2 fim AIB is the signature of the 
Yc-H band of interstellar PAHs, we could predict a band located 
at 11.0 fim if [SiPAH]^ complexes are present. Interestingly, 
such a band has been detected in many astronomical objects by 
the Short Wavelength Spe ctrometer (SWS) on board the Infrared 
Space Observatory (ISO) (H onv et al.ll200Tt iPeeters et al.ll2002l: 
van Diedenhoven et al. 

20047 We also note that an additional 
band at 6.0 fim of comparable intensity is often present. The 
presence of [SiPAH]^ complexes could account for this addi- 
tional feature: we have seen that the vc-c shift is larger for small 
[SiPAH]^ and therefore these species may be the carriers of this 
satellite AIB. Another possibility could be larger complexes con- 
taining several Si atoms. We expect the shift to increase with the 
number of Si atoms although we have not theoretically explored 
this effect yet. 

Fig. 4 shows the IR spectrum of the Red Rectangle, which 
is the brightest AIB spectrum in the sky. This region is known 
to be a peculiar C-rich nebula, where the low m etal abundances 
are expected to be linked to depletion onto dust dWaelkens et al.l 
1996). In this object, the additional band at 6.0 fim, the structure 
in the blue wing of the 6.2 fim AIB and the 11.0 fim band are 
clearly visible on the spectrum retrieved from the ISO archive. 
The strong 7.8 fim band observed in th is object i mplies that larg e 
ionized PAHs exist dBauschhcher et al . 2008; Jo blin et al.l2008i) . 
The latter can offer a large surface for multiple Si atoms adsorp- 
tion, which strengths the possible assignment of the 6.0 fim band 
to the resulting complexes. 

The region of the 6.2 fim AIB was fitted by three Lorentzian 
profiles centered at 6.0, 6.2 and 6.3 fim after a linear contin- 
uum substraction. The integrated intensities of the 6.0 and 6.2 
fim bands were found to represent 13% and 22% of the total in- 
tensity respectively. Considering that the vc-c intensity is not 
reduced by the coordination, 35% of PAHs would be involved 
in such species. Taking into account that the relative abundance 
of C in the frame work of interstellar PAHs [C/)a//]/[H] is equal 
to 6.5 X 10"^ (iJoblin et al.li 19921) and that PAHs are constituted 
by an average of 80 C atoms, that leads to a relative abundance 
[PAH1/[H1 = 8 X 10 '^ in agreement with the value given by 
Kaufman et alj (l2006l) . With regard to the cosmic abundance of 
Si ([Si]/[H] = 3.2 X 10"^ from Asplund et al. 2005), we find that 
the silicon abundance which would be involved in these species 
corresponds to an order of magnitude of 1%. 
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We compared the areas between the 11.0 and the 11.2 yum 
bands using the same method. The area of the 11.0 fim band rep- 
resents 6% of the total area composed of the 11.0 and 11.2 fj.m 
bands. This value seems to be consistent with the higher value 
extracted at 6.2 fim taking into account that (i) a redshifted coun- 
terpart which falls in the wing of the intense 1 1 .2 fim AIB is ex- 
pected at ~ 1 1 .4 yum and (ii) the coordination geometry does not 
lead to a systematic activation of the 1 1 .0 fim vibration mode. 

5. Conclusions 

iTielensI (Il998h proposed that some silicon is locked up in a dust 
component which differs from the silicates and whose bind- 
ing energy Eh is ~2 eV. Our thermodynamic results show that 
[SiPAH]^ complexes are good candidates with a minimum bind- 
ing energy Eh of 1 .5 eV. The formation of such species in the 
ISM is expected to occur by radiative association of Si^ with a 
PAH° molecule, which is an exothermic process liberating 3.0 
eV. 

The comparison of the h-iOf^m/ho-iS/jm ratio for PAH^ and 
their Si-complexes illustrates the charge effect induced by the 
coordination of silicon: this does not strongly affect the IR in- 
tensities of the spectra of PAH^ . Hence the nature of silicon en- 
ables the formation of cationic complexes in which the charge is 
mainly located on PAHs. That leads the [SiPAH]^ complexes to 
show IR spectra close to those of PAH^ in terms of intensities, 
along with characteristic features induced by the coordination of 
silicon. 

We propose that the presence of [SiPAH]^ species in astro- 
nomical environments can be deduced from the 11.0 fim satellite 
band. It can also provide an explanation for the blueshifted po- 
sition of the 6.2 yum AIB. Additionally, the 6.0 yum satellite band 
could be the signature of multiple Si atoms complexes. A firm 
assignment would require further experimental and theoretical 
work. We suggest that this identification does not require strin- 
gent constraints on the silicon abundance since typically 1% of 
the cosmic silicon would be attached to PAHs. 

A previous theoretical study on [FePAH]^ species has not 
revealed specific finger prints that could be u sed to identify 
these species in the ISM dSimon & Jobhnll200l . Our results on 
[SiPAH]^ provide the first spectroscopic evidence of the pres- 
ence of metal-PAH species in interstellar and circumstellar envi- 
ronments. 

Acknowledgements. We acknowledge the anonymous referee for his fruitful 
comments on the manuscript, and Paolo Pilleri for his help with the ISO spec- 
trum. This work was supported by French National Program 'Physique et Chimie 
du Milieu Interstellaire' which is gratefully acknowledged. We also thank the 
CALMIP supercomputing facility of Universite de Toulouse where DFT calcu- 
lations were performed. 



Hudgins, D. M., Bauschlicher, Jr., C. W., & Allamandola, L. J. 2005, ApJ, 632, 
316 

Hudgins, D. M. & Sandford, S. A. 1998, J. Phys. Chem. A, 102, 329 
Jaeger, J. B., D., R E., Jaeger, T. D., & Duncan, M. A. 2005, J. Phys. Chem. A, 
109, 2,801 

Jaeger, T, van Heijnsbergen, D., Klippenstein, S., et al. 2004, J. Am. Chem. 
Soc, 126, 10981 

Joblin, C, D'Hendecourt, L., Leger, A., & Defourneau, D. 1994, A&A, 281, 923 

Joblin, C, Leger, A., & Martin, P 1992, ApJ, 393, L79 

Joblin, C, Szczerba, R., Berne, O., & Szyszka, C. 2008, A&A, 490, 189 

Jones, A. P 2000, J. Geophys. Res., 105, 10,257 

Jones, R. O. & Gunnarsson, O. 1989, Rev. Mod. Phys., 61, 689 

Kaufman, M. J., Wolfire, M. G., & Hollenbach, D. J. 2006, ApJ, 644, 283 

Klotz, A., Marty, P, Boissel, P, et al. 1995, A&A, 304, 520 

Langhoff, S. 1996, Journal of Physical Chemistry, 100, 2819 

Leger, A. & Puget, J. L. 1984, A&A, 137, L5 

Peeters, E., Tielens, A. G. G. M., van Kerckhoven, C, et al. 2002, in 
Astronomical Society of the Pacific Conference Series, Vol. 267, Hot Star 
Workshop III: The Earliest Phases of Massive Star Birth, ed. P. Crowther, 
403-+ 

Pino, T., Dartois, E., Cao, A.-T., et al. 2008, A&A, 490, 665 
Pozniac, B. P & Dunbar, R. C. 1997, J. Am. Chem. Soc, 119, 10,439 
Simon, A. & Joblin, C. 2007, J. Phys. Chem. A, 1 1 1 (39), 9,745 
Simon, A., Joblin, C, Polfer, N., & Oomens, J. 2008, J. Phys. Chem. A, 112, 
8,551 

Srinivas, R., Hrusak, J., Sulzle, D., Bohme, D. K., & Schwarz, H. 1992, J. Am. 

Chem. Soc, 1 14, 2,802 
Szczepanski, J., Wang, H., Vala, M., et al. 2006, ApJ, 646, 666 
Tielens, A. G. G. M. 1998, ApJ, 499, 267 
Tielens, A. G. G. M. 2008, ARA&A, 46, 289 

van Diedenhoven, B., Peeters, E., Van Kerckhoven, C, et al. 2004, ApJ, 61 1, 928 
Waelkens, C, Van Winckel, H., Waters, L. B. F. M., & Bakker, E. J. 1996, A&A, 
314, L17 

Whittet, D. C. B., ed. 2003, Dust in the galactic environment, ed. D. C. B. Whittet 



References 

Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, ApJ, 290, L25 
Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Astronomical Society of 
the Pacific Conference Series, Vol. 336, Cosmic Abundances as Records of 
Stellai' Evolution and Nucleosynthesis, ed. T. G. Barnes, 111 & F. N. Bash, 
25-+ 

Bauschlicher, Jr., C. W. 2002, ApJ, 564, 782 

BauschHcher, Jr., C. W., Peeters, E., & Allamandola, L. J. 2008, ApJ, 678, 316 
Bohme, D. K., Wlodek, S., & Wincel, H. 1989, ApJ, 342, L91 
Bohme, D. K., Wlodek, S., & Wincel, H. 1991, J. Am. Chem. Soc, 113, 2,802 
Dunbai-, R. 2002, J. Phys. Chem. A, 106, 9809 

Dunbar, R. C, Uechi, & Asamoto. 1994, J. Am. Chem. Soc, 116, 2,466 
Frisch, M. J., Trucks, G. W., Schlegel, H. B., et al. 2003, Gaussian 03, Revision 

C.02, Gaussian, Inc., Wallingford, CT, 2004 
Hony, S., Van Kerckhoven, C, Peeters, E., et al. 2001, A&A, 370, 1030 



